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Abstract 

In the context of the minimal supersymmetric seesaw model, the CP-violating 
neutrino Yukawa couplings might induce an electron EDM. The same interac- 
tions may also be responsible for the generation of the observed baryon asym- 
metry of the Universe via leptogenesis. We identify in a model-independent 
way those patterns within the seesaw models which predict an electron EDM 
at a level probed by planned laboratory experiments and show that negative 
searches on r — > 67 decay may provide the strongest upper bound on the elec- 
tron EDM. We also conclude that a possible future detection of the electron 
EDM is incompatible with thermal leptogenesis, even when flavour effects are 
accounted for. 
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1 Introduction 



New physics has to be invoked in order to explain why neutrinos are massive and mix 
among each other, and also why they turn out to be much lighter than the other known 
fermions. Among the proposed explanations, the idea that neutrino mass suppression is 
due to the decoupling of heavy states at a very high energy has prevailed in the last three 
decades. This principle is behind the formulation of the well-known seesaw mechanism [1] . 

Besides many experiments aiming to characterize neutrino oscillations, several other 
present and future experiments are planned to search for alternative signals of lepton 
flavour violation (LFV), e.g. in radiative charged lepton decays. It is therefore of great 
importance to explore theoretical scenarios where LFV is enhanced to levels at reach of 
planned experiments. One of the most appealing theoretical frameworks where this may 
actually occur relies on supersymmetry (SUSY). Within the minimal supersymmetric ex- 
tension of the standard model (MSSM), LFV can be communicated to the SUSY- breaking 
sector by the same interactions which participate in the seesaw mechanism producing a 
misalignment between leptons and sleptons [2]. This turns out to be particularly inter- 
esting in the case where the dynamics responsible for breaking SUSY is flavour-blind, like 
in minimal-sup ergravity scenarios (mSUGRA). 

The phenomenological aspects related to supersymmetric seesaw mechanisms have 
been widely studied in the literature, especially those which concern LFV decays [3, 4, 
5, 6, 7]. In particular, most of the analysis have been based on the extension of the 
MSSM particle content where right-handed (RH) heavy neutrinos (responsible for the 
neutrino mass suppression) are added. RH neutrinos may play a role in the generation 
of the baryon asymmetry. If their couplings to the ordinary charged leptons are complex 
and therefore violate CP. then their oTit-of-cquilibrium decays may induce the observed 
baryon asymmetry of the Universe, = (0.87 ± 0.03) x 10~^° [8], via the leptogenesis 
mechanism [9]. In the simplest thermal scenario, the lightest heavy Majorana neutrino is 
produced after inflation by thermal scatterings and subsequently decays out-of-equilibrium 
violating CP and lepton number. Flavour effects play also a role in thermal leptogenesis 
since its dynamics depends on which charged lepton Yukawa interactions are in thermal 
equilibrium and, therefore, the computation of Yb depends on which temperature regime 
the heavy Majorana neutrino decays occur. We will come back to this point later. 

Besides flavour violation, also CP-violating (CPV) effects may be induced in the soft 
SUSY breaking terms by the CPV seesaw interactions. These may lead to relevant con- 
tributions to the charged lepton electric dipole moments (EDM), in particular that of 
the electron [10, 11, 12, 13, 14]. The present upper bound on the electron EDM is 
de < 1.6 X 10~^^ e cm at 90% C.L. [15]. Within three years, the Yale group plans to reach 
a sensitivity of about 10^^^ e cm [16] and hopefully go down to 10^^^ e cm within flve 
years. A more ambitious proposal, based on solid-state physics methods, exists to probe 
de down to 10~^^ e cm [17]. 
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The above discussion shows that, by extending the MSSM particle content with heavy 
neutrino singlets, a new window is widely opened into the investigation of new effects 
in flavour physics and cosmology. Ultimately, one expects to be able to relate various 
phenomena like neutrino oscillations, LFV, leptogenesis and EDMs. This programme 
is of extreme importance and may shed some light over the origin of neutrino masses 
and mixing and its relation with possible new physics observations. In what follows, 
we consider the scenario where the (dominant) LFV and CPV effects in the SUSY soft 
breaking sector are exclusively generated by the seesaw Yukawa interactions. Under these 
assumptions, we identify the Yukawa structure leading to an observable de and at the 
same time compatible with the present bounds on LFV decays and show that there is an 
upper bound on de coming from negative searches of the r — > 67 decay. This indirect 
bound can be even stronger than the direct-search limit. We will also study the impact 
on thermal flavoured-leptogenesis for cases where de is at hand of future experiments, 
showing that a positive detection of de is incompatible with thermal leptogenesis. 

The paper is organized as follows. In Section 2 we review the issue of EDM and 
LFV within supersymmetric seesaw models and discuss the bound on de provided by the 
T — > 67. In Section 3 we give a short review of thermal leptogenesis, including flavour 
effects. In Section 4 we discuss the patterns of neutrino Yukawas giving rise to a large 
de, but still compatible with present bounds on — >• ej, r — >• 67 and r —>■ /I'y. The 
conclusions are drawn in Section 5. 

2 EDM and LFV from seesaw Yukawa couplings 

We consider the framework of the minimal supersymmetric standard model (MSSM) 
extended with three heavy Majorana neutrino singlets N. Working in the basis where the 
charged lepton mass matrix is real and diagonal, the seesaw [1] interactions at high 
energy are described by the following superpotential terms 

Wss^NY,LH^ + ^NMN , M = diag(Mi, M2, M3) , (1) 

where Y^, is the Dirac neutrino Yukawa coupling matrix and Mj stand for the (positive) 
heavy Majorana neutrino masses ordered as Mi < M2 < M3. As usual, L denotes the 
lepton doublets and the hypercharge 1/2 Higgs superfleld. Integrating out the heavy 
neutrino singlets, after electroweak symmetry breaking one obtains the effective Majorana 
neutrino mass matrix 

jji^ff = —Yjyv't = U*mU^ , m = diag(mi, m2, ma) , (2) 
M 

where U is the leptonic mixing matrix (parameterized in the usual way, see Eq. (22)), rrii 
are the (positive) effective neutrino masses, — v sin /5 with v — 174 GeV and tan (3 — 
{Hu)/{H,). 
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Because of RGE running from high to low energy scales, the seesaw Yukawa couplings 
potentially induce flavour and CP violations in slepton masses. In the following we review 
their impact by assuming universal and real boundary conditions at the Planck scale Mpi: 
rrix — rnl for the soft scalar masses, Ax — aoYx for the trihnear A-terms and Mi — M1/2 
for the gaugino mass parameters. In addition, we assume that the //-term is real. Under 
these circumstances, the flavour and CP violations in the low-energy slepton masses arise 
from those present in the seesaw Yukawa couplings through RGE running. 

The most relevant flavour misalignment is induced by the seesaw Yukawas in the slep- 
ton doublet mass matrix m| [2] . Following the widely used mass insertion approximation, 
the flavour- violating entries of m| can be parametrized by 

~ - (47r)2 ml ^ ' 

where is an average doublet slepton mass and Cjj encodes the dependence of on 
the Yukawa interactions: 

^^. -E^S ' = (4) 
k ^ 

For later convenience, we have isolated the contribution of each iVfc to Cij denoting it by 
Cf-. Notice that in Eq. (4) has to be evaluated at Mpi. The induced LEV decays are 



BR(£, ^ £,7) = BR(£, ^ tjVjV,) Fb tan^ (3 , (5) 

where Fb is an adimensional function of supersymmetric masses which includes the con- 
tributions from chargino and neutralino exchange - see e.g. [6] and references therein. 

From the above considerations, it is clear that the present experimental limits on LEV 
decays can be translated into upper bounds on \Cij\, denoted by CfJ^, which depend on 

tan/? and supersymmetric masses [5]. In mSUGRA, once the relation between qq and 
mo has been assigned and radiative electroweak breaking required, the supersymmetric 
spectrum can be expressed in terms of two masses, e.g. the bino mass Mi and the 
average singlet charged slepton mass rfiR at low energy (we recall that Mi ^ O.AM1/2 
and m|. ^ + {i.l^M^^. For definiteness, we will take from now on = + Mi/2- 
Considering in particular the point P = {Mi^rfiR) = (200, 500) GeV, for which the SUSY 
contribution to the anomalous magnetic moment of the muon SUf^ is within the observed 
discrepancy between the SM and experimental result for tan/3 > 35, we obtain 

C^^^BxlO-'-^ , Cgtc^O.S^ , ^3^=:^^ . (6) 
tan/3 tan/3 tan/3 ^ ^ 

The strongest constraint comes from fi — 67, although also those from r — 67 and 
T ^ /I'j become non-trivial for tan /3 > 10, as they imply that couplings are at most 
of 0{1). This, in turn, supports the use of the perturbative approach. In Fig. 6 of the 
Appendix we display the dependence of Cfj' in the plane (Mi, fhu). 
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There are two potentially important sources of CP violation induced by the seesaw 
Yukawa couplings in the doublet slepton mass matrix. One is associated to the flavour- 
conserving (FC) 74-terms [10, 11, 12], while the other, generically dominant for tan /? > 10, 
is mediated by flavour- violating (FV) 5's [12, 13]. The corresponding contributions to 
lepton EDMs are given by 

^ [e cm] = Fd mi, Im(ai) , ^ [e cm] = F'^ ^ tan 13 Im(5^^m^ 5^^)u , (7) 

where F^i, F^ (with dimension of mass~^) are functions of the slepton, chargino and 
neutrahno masses - see e.g. [6] and references therein. These sources of CPV are 

Tr^^n ^ ^"Q tFC Jr^fxRRm sLL\ Sm^ (6mg + 2ag)(6mg + 3ag) ml tan^ (3 py 

(8) 

where 

fc>fc' k>k' 

Adopting a notation which renders more explicit the link with LFV decays and deflning 
In^ = ln(M„/M6) for short, one has [14] 

7^')^° = ^ MC'^O,, , it'^^-" ^Inl I J C'^^^C') , (10) 

l^fe' V "^r J a 

with ln2 = In^, In^ = ln?(l - 2 In^ In? / In^ In^') , In^ = ln?(l - 2 In^Hn^ / In^Hnf ) . Notice 
that these contributions arise as an effect of a splitted spectrum of right-handed neutri- 
nos and would vanish in the case of a degenerate spectrum [10, 18]. It turns out that 
the seesaw-induced contributions to and d^. are below the planned experimental sen- 
sitivities; on the contrary the seesaw-induced contribution to might be at the level of 
planned experiments. The present experimental upper limit, d'^f^ = 1.4 x 10^^'' e cm, can 
correspondingly be translated into upper bounds on and with a dependence 

on tan/3 and supersymmetric masses - see [12, 14] for more details. Within the mSugra 
framework, we display these bounds in Fig. 6 of the Appendix. In particular, for the point 
P introduced before, the upper bounds are 

^ub,FC ^ ^ub,FV ^ ^q2 ( ' 

\tanPj ^ ' 

For comparison, it is useful to estimate the upper allowed values of Jf and Jj^ assuming 
perturbativity: by allowing the relevant Yukawa couplings of Yy to be of OiV) and the 
logarithms to be large enough, from Eq. (10) one obtains Jf" < 50 and Jj^ < 300. 
This means that in point P the FC-type seesaw-induced contribution, dg^, is below the 
level of the present experimental sensitivity; an improvement by two orders of magnitude, 
pushing dg^*' at the level of 10~^^ e cm, would be needed to start testing it. On the 
contrary, the FV-type seesaw- induced contribution, d^^, might have already exceeded 
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the present experimental bound in P for tan/? > 35 (and for even smaller tan/3 if the 
supersymmetric masses are taken to be smaller than those in P). 

However, as we now turn to discuss, also the experimental limit on r — > 67 provides 
an upper bound to the seesaw- induced Ij^ and Ij^. At present, such indirect bounds are 
even stronger than those from direct searches of the electron EDM. Once the constraint 
from T — > 67 is taken into account, it turns out that both the FC and FV-type seesaw 
contributions have to be sizeably smaller than what would be allowed assuming pertur- 
bativity. The argument is based on the fact that, barring cancelations, \C^j\ < C^J^ for 
each k. In Eq. (10), the strong bound from // — 67 makes the terms involving to be 
negligibly small. If large, /J*" and /J^ are then proportional to the same combination of 
Yukawas [14]: 



In the discussion of the next sections, it will turn out that the only relevant contribution 
comes from {k,k') = (3,2). Therefore, in what follows, we will focus on /i^^-*^*^ ln2 ^ fa 
^(32)FV ^ \ullm{Cf^Cf^). Requiring lm{C^^C^'^*) < C^g^, one thus obtains an indirect 
upper bound from r — > 67, to be denoted by /J. Considering in particular the point P, 
the indirect upper bounds from r — > 67 are 



As a comparison, allowing the relevant Yukawa couplings Y,^ to be (at most) of 0{1) 
and the logarithms involved to be large, one would have lm{Cf^CfS^) < 50. This means 
that, for tan/? > 5, both /J^*" and I^^^ are smaller than what allowed by perturbativity. 
Typically, 10 JJ^^ - If^ - 500 for tan/3 = 5, while Wlf^ - If^ - 5 for tan/3 = 50. 
The indirect bounds from r — > 67 are thus stronger than those from direct searches of 
de, shown in Eq. (11). Taking P and representative sphttings between M3 and M2, we 
display in Fig. 2 the ratio dl/dfP = /jFC/jub.FC ^ jrFV/jub,FV ^ function of tan/3. It 
turns out that the indirect upper bound from r — > 67, dl, is stronger than the present 
direct one, dl^^, by about two orders of magnitude"*^. Planned experiments lowering the 
bounds on d^ by about three orders of magnitude will then provide sensible tests of the 
seesaw-induced effects. Clearly they will test the seesaw models that maximize Jf ^ (hence 



also Jf^), I.e. the models where: 1) M3 > 10 M2; 2) lm{Cf^Cf^) ~ C^f. The latter of 



course also implies that r — 67 is close to the experimental bound. 

^As a consequence, a potential discovery of the electron EDM within an order of magnitude from the 
present sensitivity should not to be interpreted as due to the seesaw-induced effects - this of course holds 
for mSugra and point P. 




(12) 




(13) 
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Figure 1: Ratios min(/J^c^ 50/ Inf)/!^'^^ (dashed) and mm{If^, In^ (dash- 
dotted) evaluated in P as a function of tan/3. Their sum is d'l./d^f^^ (sohd). We take 
M2 = 5 X 10^^ GeV and M3/M2 = 10^ 10^, 10 from top to bottom, respectively. 



3 Thermal leptogenesis including flavours 

In this section we provide a short review of leptogenesis which will be useful to establish the 
link between the high energy CP-violation responsible for the generation of the baryon 
asymmetry through thermal leptogenesis and the one responsible for the EDM of the 
electron. 

Thermal leptogenesis [19, 20, 21] takes place through the decay of the lightest of 
the heavy Majorana neutrinos if these are present in the early Universe. The out-of- 
equilibrium decays occur violating lepton number and CP, thus satisfying the Sakharov's 
conditions [22]. In grand unified theories (GUT) the heavy Majorana neutrino masses are 
typically smaller than the scale of unification of the electroweak and strong interactions, 
Mgut — 2 X 10^^ GeV, by a few to several orders of magnitude. This range coincides 
with the range of values of the heavy Majorana neutrino masses required for a successful 
thermal leptogenesis. 

We will account for the flavour effects which appear for Mi < 10^^ GeV and have been 
recently investigated in detail in [23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34] including the 
quantum oscillations/correlations of the asymmetries in lepton flavour space [25, 34]. The 
Boltzmann equations describing the asymmetries in flavour space have additional terms 
which can significantly affect the result for the final baryon asymmetry. The ultimate 
reason is that realistic leptogenesis is a dynamical process, involving the production and 
destruction of the heavy RH neutrinos, and of a lepton asymmetry that is distributed 
among distinguishable lepton flavours. Contrarily to what is generically assumed in the 
one-single flavour approximation, the AL = 1 inverse decay processes which wash-out the 
net lepton number are flavour dependent, that is the lepton asymmetry carried by, say, 
electrons can be washed out only by the inverse decays involving the electron flavour. The 
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asymmetries in each lepton flavour, are therefore washed out differently, and will appear 
with different weights in the final formula for the baryon asymmetry. This is physically 
inequivalent to the treatment of wash-out in the one-flavour approximation, where the 
flavours are taken indistinguishable, thus obtaining the unphysical result that, e.g., an 
asymmetry stored in the electron lepton charge may be washed out by inverse decays 
involving the muon or the tau charges. 

When flavour effects are accounted for, the flnal value of the baryon asymmetry is the 
sum of three contributions. Each term is given by the CP asymmetry in a given lepton 
flavour properly weighted by a wash-out factor induced by the same lepton number 
violating processes. The wash-out factors are also flavour dependent. 

Of course, since we are dealing with the MSSM, we have to consider the presence of the 
supersymmetric partners of the RH heavy neutrinos, the sneutrinos Ni {i = 1, 2, 3), which 
also give a contribution to the flavour asymmetries, and of the supersymmetric partners 
of the lepton doublets, the slepton doublets. Since the effects of supersymmetry breaking 
may be safely neglected, the ffavour CP asymmetries in the MSSM are twice those in 
the SM and double is also the possible channels by which a lepton ffavour asymmetry 
is reproduced. However, the AL = 1 scatterings washing out the asymmetries are also 
doubled and the number of relativistic degrees of freedom is almost twice the one for 
the SM case. As a result, introducing new degrees of freedom and interactions does not 
appreciably change the flavour asymmetries. There are however two other and important 
differences with respect to the nonsupersymmetric thermal leptogenesis SM case. First, in 
the MSSM, the flavour- independent formulae can only be applied for temperatures larger 
than (1-l-tan^ /3) x 10^^ GeV. Indeed, the squared charged lepton Yukawa couphngs in the 
MSSM are multiplied by (1-l-tan^ Consequently, the ji and r lepton Yukawa couplings 
are in thermal equilibrium for (1 + tan^ /3) x 10^ GeV < T < (1 + tan^ /3) x 10^ GeV and 
all flavours in the Boltzmann equations are to be treated separately. For (1 + iw? j3) x 
10^ GeV ^ T <^ (1 -|- tan^/9) x 10^^ GeV, only the r Yukawa coupling is in equilibrium 
and only the r flavour is treated separately in the Boltzmann equations, while the e and 
11 flavours are indistinguishable. Secondly, the relation between the baryon asymmetry 
Yb and the lepton flavour asymmetries has to be modifled to account for the presence of 
two Higgs fields. Between (1 -Ftan^/3) x 10^ and (1 -Ftan^/3) x 10^ GeV, the relation is 
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93 ^ \ /19^ \ /19 



' ^ 110 y ^ ' V30 ^/ ' V30 



(14) 



where the fiavour lepton asymmetries are computed including leptons and sleptons. Be- 
tween (1 + tan^ j3) X 10^ and (1 tan^ /3) x 10^^ GeV, the relation is 



Sit/* 



541 _ \ /494^ 
'^61 7 ' V761 



(15) 



where the number of relativistic degrees of freedom is counted by — 228.75. The 
observed value is Yb = (8.7 ± 0.3) x 10~^^ [8]. Let us explain the various terms entering 
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in Eqs. (14) and (15). The CP-asymmetry in each flavour is given by [23, 26, 25, 

3Mi ^"^{^pp^T^Tui^UipRipRi^ 



S/3 "^/3 |-Rl/3| 

where the (in general complex) orthogonal matrix R [4] is defined as 



= RVAUK 



(16) 



(17) 



Similarly, one defines a "wash-out mass parameter" for each fiavour £ [23, 26, 25, 

2 



2 sin^/3 X 10-3 eV 



//(Ml) 



, me = 



Ml 



(18) 



The quantity me parametrizes the decay rate of A^i to the leptons of fiavour i. Further- 
more, in Eq. (15), €2 — ee + e^^ and m2 = rhe + mp^. Finally, 



V {me) ~ 



me 



8.25 X 10-3 eV 



0.2 X 10-3 eV 



me 



-1.16 



(19) 



parametrizes the wash-out suppression of the asymmetry due to AL = 1 inverse decays 
and scatterings. Notice that the wash-out masses TO2 and rhr in Eq. (14) are multiphed by 
some numerical coefficients which account for the dynamics involving the lepton doublet 
asymmetries and the asymmetries stored in the charges = (1/3)5 — [28]. 



Maximizing with the constraint from /x 
and impact on thermal leptogenesis 
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In this section our goal is to identify those patterns within the seesaw models which 
predict de at hand of future experiments, while keeping the prediction for LFV decays, in 
particular // — > 67, below the present bound. In this respect, the analysis we perform is 
model-independent. From the discussion of the previous section, it is clear that to have de 
at hand of future experiments, the quantities 1^ defined in Eq. (9) have to be maximized. 
Our aim is also to understand whether these models can explain the observed baryon 
asymmetry of the universe via Ni decays. Our work differs from that of Ref. [11] in two 
respects. First, our analysis includes also the FV contribution to the electron EDM, which 
is dominant for large values of tan (3 and, secondly, we include fiavour effects in computing 
the baryon asymmetry from leptogenesis. 
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We choose to work with not too heavy RH neutrino Ni, say Mi in the range (10"^° — 
10^^) GeV. This choice has the advantage that, barring accidental cancelations, the cou- 
plings of A^i are rather small, Y^ij < a; = C(10~^). This has important consequences: i) 
C21 does not exceed the — > 67 bound; ii) the eventual contributions to dg from 7^^^^'^^^^ 
(both FV and FC) are much smaller than the limit inferred from experiment and can be 
neglected in first approximation. The only potential source of a large is therefore 

7f In^ ^ If^^^ ^ h \nf \nl , h = MY:,,Y^,,Y^,,Y:,,) , (20) 

which is maximized requiring that: 

1) the splitting between M3 and M2 is large; 

2) the four Yukawas appearing in Eq. (20) are of 0{1), but satisfy the bound from r — > 67. 

Barring conspiracies between Cf^ and Cf^, to suppress /i — > 67 below the present 
bound one must require both to be smaller than C21. Since oc ^jy*fc2^i/fcu order to 
keep li^^-* large, one has to impose |5^i/32|, l^^i/ssl ^ 2/ ~ / 5. Explicitly, we are looking 
for textures of the form 

/ <x \ 
n= Oil) <y Oil) , (21) 
\Oil) <y Oil)) 

where, as already mentioned, x = (9(10~^) ensures that C21 satisfies the upper bound on 
II — > ej. Notice also that this pattern enhances r ^ ej while suppressing r fi'j. 

It is convenient to exploit the complex orthogonal matrix R [4] introduced in Eq. (17) 
and adopt a standard parameterization for the MNS mixing matrix 

U = 023(^23) 013(^13) r| 012(^12) X diag(e-l/^ e-^/^ 1) , (22) 

where Ts = diag(l, 1, e**^), Oij = [(Qj, %)(— Sjj, Q^)] with Cij = cos% and Sij = sinOij. 
The Dirac and Majorana type phases were denoted by 6 and ai 2, respectively. We adopt 
for R a parameterization in terms of three complex angles 9[j: 

(^R ^R ^R qR qR qR qR I ^R ^R ^R 
''12'-13 "'l2'-23 '-12*23*13 *12"'23 ~^ '-12'-23"'l3 
-44 cf24 + 5f244 44-444 
r^R r^R „R „R „R 
•^13 '^23*^13 *-23*^13 

(23) 

As already discussed, the 67 constraint is implemented by requiring 1 1^1,32 1, |y^22| % 
y. In the approximation y — (which turns out to be very satisfactory), this means that 
R must satisfy respectively the conditions 

Ri3V^u;^ + Ri2V^u;2 + RiiV^iu;^ = o , i = 3,2 . (24) 

Simultaneously we want to enhance 

M3M2, 



-Im 



iRVrhU^)l^iRVrfiU^)33iRV^U^)2iiRVrfiU^)*2^ . (25) 
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A potentially interesting contribution to the electron EDM requires M3 and M2 to be 
as large as possible and, as already stressed, enough splitted to ensure a large ln2. This 
also explains why we cannot rely here on resonant leptogenesis [35]. One could have, for 
instance. Mi ~ M2 <^ M3. But since Mi has to be small (for ji — > 67), this would imply 
small M2 and consequently a suppressed 1^. Instead, the pattern Mi <S M2 ~ M3 implies 
lug -C 1. Moreover, any baryon asymmetry resulting from the resonant decays of N2 and 
iVs can be washed out by A^i (see however [36] and especially [37]). Finally, the situation 
Ml ~ M2 ~ M3 would lead to a very suppressed [10, 18]. 

As we are going to discuss in the next section, the structure of in Eq. (21) deter- 
mines, to some extent, the form of R. From a qualitative point of view, one can already 
put forward some guesses. The requirement of a large atmospheric angle implies that for 
some Nj it should be Yi,j2 ~ ^i/js, but from Eq. (21), this can happen only for j = I. In the 
case of hierarchical light neutrinos, this implies dominance of ms by Mi, while M3 and M2 
will be mainly associated to the lighter masses. Since we aim at M2 << M3, this means 
that M2 and M3 will dominate respectively m2 and mi. According to the interpretation 
of i? as a dominance matrix [38] , these considerations in particular determine the first row 
of R, which is relevant for leptogenesis, to be approximately (0, ^/m2/rr^3, 1 — 0(777.2/7713)). 

In the following we treat separately the cases of normal and inverted hierarchy for 
light neutrinos. 



4.1 Normal hierarchy for Ught neutrinos 

For normal hierarchy (NH) we take ~ m@ = \/ml — ml, 1712 ~ 7770 = ^/ ml — ml and 
leave 777i <^ 7772 undetermined. It is convenient to rewrite the conditions (24) assuming 

Ri3 ^ 0, 

-R,2 - Ri\ ~ 

—^ = t+-^t, (26) 
where the quantities t and t depend only on the neutrino parameters at low energy. 



1713 U23 _ frn^ 2 ^23 



7772 C^22 V ^2 C12 



(27) 



mi Uoi Pmi 



7772 1^22 V 

with tij = tan 9ij. In particular, for a light neutrino spectrum with NH and taking into 
account the present neutrino oscillation data, \t\ ~ 3 while \t\ <^ 1. 

Using the above expressions and introducing the reference scale M@ = v'^/m@ ~ 
5 X 10^^ GeV for the right-handed neutrino masses, the Yukawa couplings of the second 
and third rows of are given by 



1^2. = ^/|^(i^23u:+i?2r.pw: 

M@ V V "^3 
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^--V1H'+^"^*')' '''' 

with 

= C/i3 - , Wi = [/a-^C/i2. (29) 

<^22 t722 

Notice that U2 and W2 correctly vanish, as required by the conditions 1^1,22, Yu32 — 0. We 
also anticipate that, for mi = 0, one has I 4 oc Im(|Uip|U3p) = which, in the present 
framework, leads to a very suppressed de- For the sake of the following discussion, we 
report the approximate expressions for Ui,3 and Wi,3 at first order in S13: 

Ui = -t23tl2(l + ^23^125136*') + 5136"'^ , Wi = 6^^ — {1 + 5136^3^12) , 

C12 

U3 = — (l + si3e''i23ii2) , W3 = -si3e^('+^)^— . (30) 

C23 C23C12 

As an example, for tri-bimaximal mixing one has: U3 = V2, Ui = -1/v^, W3 = and 
Wi = 73/2 e^"'^^- 

Prom the above expressions, and taking into account the smallness of W3 = 0(C/e3), 
it turns out that large yj,23 requires i?23 = C^(l) and M2 = 0(M@). Then, the condition 
of large splitting between M2 and M3 implies M3 >> M@ while large 1^1,33 requires -R33 = 
0{^M@/Ms) < 1. Exphcitly, i?33 ~ ^23^13 t»ut, since the condition for i — 3 in Eq. (26), 
^23 — ^ ~ ^13^*^23) naturally suggests c^^ to be largc^, the parameter to be suppressed is 
rather c^^. We then define cfg = % < 1 and expand at first order in x- In particular, 
from Eq.(28), we have c^^x ~ V M@ /M^Yi,^^ /ILJ3 . The conditions in Eq. (26) get now 
simplified: 

y^R ^.R I r,R. r>R. 1 

g^V^«c« + ° (^^) = * = + -^t. (31) 

*"12*23 "I" *12*"23 /V *-23 

One can then envisage two relevant cases according to the value of 6^2 ■ 

• If sfs = e < 1, Eq. (31) gives t^^ -1/f + e and ff^ {t + 1/t) XC23. The latter can 
be rewritten as 

V m3 V M@ S12 

Here, mi cannot be arbitrarily small, otherwise M3 would exceed the GUT or Planck 
scale. At first order in e one has Ri —c + es, s + ec, so that 

X c s 

i? = I 7s -c\ +0(e^ex,x'). (33) 
-1 xc XS 



^Clearly, c|| is small only when ^23 ~ 7''/2, but in this case the condition for i = 3 in Eq. (26) forces 
cfg to be still small, cfl ~ —t. 
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As expected, since c ~ 1/i < 1 and s ~ 1 — 1/(2 P), this structure means that 
ma ~ m@ is dominated by the hghtest right-handed neutrino A^i, while m2 ~ 
by N2. Instead, mi is dominated by the heaviest, A^3, which decouples the more x 
is small. Notice that such a precise determination of the first row of R allows to 
predict leptogenesis, as we are going to discuss. The Yukawa couplings relevant for 
dp are 



Yu3i — 



Ms 



sx^-. r^W*) , Y,2i 

V "^3 



M2 



-cU* + C»(x)W — W; 



(34) 



so that 



M2 M3 
M@ M@ 



mi 



Im -sxWiUtlUal" + sxWaU^lUil' + 



mi 
ms 



— 1^33^1^31^21^231 



^ sm a2 + ^(513 sm d) 



(35) 

and displays a suppression 



The first term is the only one present in the limit | Ues 
by a factor m^/m® ~ 0.17. For large values of |f/e3| and 8 ~ ±7r/2 (large Dirac-type 
CP violation at low-energies), the second and third terms can be dominant. 

If = e < 1, Eq. (31) gives ifg i (1 - ) and Fsi (i + l/t)xc2^- The latter 
can be rewritten as 




ma V M@ 



\et 



S23| 
Sl2 



i/33 



0(6). 



(36) 



In this case mi is allowed to be small and the limit mi — > can be applied. At first 
order in e one has ~ s — ec, ~ c -|- es, leading to 



R 




0{e\ex.x'). 



(37) 



Again, dominance of A^i and N2 has been obtained respectively for m® and m©, 
while mi is associated to iVs, which decouples the more x is small. The Yukawa 
couphngs are now 



Ms 



m3 



M2 

M® 



-c-u*+o(x)J^w: 

m3 



(3^ 



and, consequently. 



c 




= £ I ^i/33^i/31 ^1/21 ^1/23 1 



Im -cxWiU^lUal^ + cxWaU^jUil^ + 



mi 
ms 



^23 . ^2 , . rx 

-— sm -— -|- 0{Si2, sm 0) 
C12 2 



(39) 
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which is suppressed the more e is small. Indeed, as already mentioned, I4 vanishes 
in the limit mi = (namely e = 0), as a consequence of the alignment between the 
second and third rows of Y^,, which ensures a vanishing eigenvalue. The remaining 
contributions from J^^^^'^^'^^ do not vanish but, as stressed before, are negligible even 
considering future sensitivities. 

Summarising, it turns out that a large d^. is obtained only under the conditions of the 
first case, namely if M2 ~ M@ and the relation mi/m@ ~ M@/M^ holds. The contribution 
induced by 5 could even exceed the one associated to 0:2; the phase ai plays no role. Notice 
in particular that values of mi < 0.5 x 10~^m3 are incompatible with an experimentally 
relevant de- 

We remark that in all the expressions above, the Yukawa couphngs, the elements of 
the matrix U and the light neutrino mass eigenstates are consistently evaluated at Mpi. 
For instance, in the case of NH, the effect of running is such that m is rescaled by a 
numerical factor, while the parameters of U do not change significantly. We illustrate this 
point by considering a numerical example for the first case discussed previously. If we 
take the low energy neutrino spectrum consistent with the present-day neutrino oscillation 
data, m = (0.96 x lO^^ 0.17, l)m@, then, for tan/5 = 30, (M3,M2) = (102,2)M@, Mi = 

3 X 10^"*^ GcV, a2 = vr/2, 6 = 7r/2 and any «i, at the Planckian scale we obtain^ m = 
1.8(10^^,0.166, l)m@. Correspondingly, the mixing angles change from 623 = 47°, 6*12 = 
35.2°, ^13 = 9.7° to 023 = 45°, ^12 = 35°, ^13 = 10°. With these choices one obtains: 
7jv = 10.2, = 1.35, C31 = 1.5, C32/8 = C21 = 10-^ These values are slightly changed 
taking different values of tan /3, which enters only through the running effects. This means 
that, in the point P and for any tan (3, the present example gives d^^ = 1.3 x 10~^dl^^, 
which may escape detection. Since d^^ and LFV decays explicitly depend on tan f3, we 
need to specify it. Taking for instance tan (3 = 30, one has d^^ = 0.022 (if P, which is at 
hand of future experiments. With tan/9 = 30, the BRs of r ^ 67, r ^ and /i — 67 
turn out to be smaller than their corresponding experimental limits by factors of 0.8, 

4 X 10~^ and 0.02, respectively. 

The left-hand side of Fig. 2 shows the dependence of /f ^ and Yg on a2 for a given set 
of values for S, keeping for the remaining parameters of U and light and heavy neutrino 
masses the same set as before - in particular we recall that Mi = 3 x 10^^ GeV and 
tan/3 = 30, but the dependence of y,^(Mpi) on tan/3 is mild. On the right-hand side of 
the same figure, we show the behaviour of /f^ and Yb in the {a2,S)-p\ane. Hence, /f^ 
displayed in Fig. 2 slightly changes considering other values of tan/3. The same applies 
separately for each asymmetry e^, but not to Yb, for which two different regimes can be 
identified according to the value of tan /3. Indeed, having Mi = 3 x 10^^ GeV, in the 
case that tan/3 > 10 the baryon asymmetry is generated in the range of temperatures 

■^In practice we do the opposite, namely we assign U and rescale m by a numerical factor at high 
energy and check that, when evolved at low energy, the parameters of U and the neutrino spectrum are 
still within the experimental window. 
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Figure 2: Dependence of /f^ (left upper-plot) and Yb (left lower-plot) on 02 for a given 
set of values for 5 and tan/5 > 10 (/i and r in equilibrium): we take 5 = 0, 7r/2, vr, 37r/2. 
On the right-hand side we show the behaviour of /J^ (upper plot) and Yb (lower-plot) in 
the (0^2, (5)-plane. For the choice of the remaining parameters, see the text. 

for which all the lepton flavours, but the electron one, are in thermal equilibrium and 
Eq. (14) applies, see Fig. 2; in the case that tan/5 < 10 thermal leptogenesis takes place 
in a range of temperatures where only the tau flavour is in thermal equilibrium and, in this 
case, one has to apply Eq. (15), see Fig. 3. From the comparison of the upper and lower 
plots of Fig. 2 one can see that there are regions in the (^2, 5) parameter space where the 
seesaw- induced effects can lead to values of de within future experimental sensitivities and 
Yb is compatible with its observed value. For tan /5 < 10 the baryon asymmetry is on the 
contrary below its observed value (the horizontal grey line) and one would need slightly 
larger values of the mass Mi to obtain a value of the baryon asymmetry consistent with 
observation. This seems to be a generic conclusion. Also notice that, in this temperature 
regime, Yb does not depend strongly on 5. 

Fig. 5 shows the dependence of /f^ and Yb on nii/m^ (left and right panel respec- 
tively), keeping mi/m^ = M@/M^ as in the previous example, as well as the same set 
for the remaining parameters. We also display (dashed line) the baryon asymmetry com- 
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3nl2 



to 71 




Figure 3: The same as in Fig. 2 for taking tan/3 < 10 (only r in equilibrium). 

puted in the one-single flavour approximation. We see that there is almost one order of 
magnitude difference. 

4.2 Inverted hierarchy for Ught neutrinos 



For inverted hierarchy (IH), we take m2,i ~ m@ and leave undetermined. It is con- 
venient to redeflne the matrix R according to: R = O12O13O23 x antidiag(l, 1, 1). This 
simple exchange of the flrst and third columns of R allows to carry on the discussion in 
a parallel way to what done before. Indeed, assuming also Rn 7^ 0, the conditions (24) 
become 

Ri2 



Rii 



~ Ri% - 



(40) 



where t and t are deflned as in Eq. (27), but now we have \t\ ^ 1/v^ while \t\ <^ 1. 

Using these expressions, the Yukawa couplings of the second and third rows are given 

by 



u2i 



Mo 



M, 



@ 



R' 



21' 



+ R 



M, 



R 



31' 



+ R 




(41) 



which, as before, correctly vanish for i = 
that to maximize de we require M3 ^ 



2. Large Fy2i requires M2 = (9(M@), so 
. Instead, large F^si is obtained provided 



Ri3 = 0(^M@/M3) < 1. Exphcitly, R31 = cfgcfg, but the condition for i = 3 of Eq. (40), 



''23 



t — ^13^/0^3, suggests that is quite large. We rather suppress cfg = x < 1 



and expand at flrst order in x- Iii particular, now 0^3 x 
conditions in (40) get simplifled: 



^.R _i_ r^R ^R 
'-12'-23 "T ■='12'='23 



M C.R 



.R ^R 



"'^12 -^23 + '^12'^23 



+ oixt) = t 



4 + 



y/M@/MsY^si/Wl and the 



(42) 



XC23 
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It is useful to distinguish between two relevant cases. 



• If sfs = e < 1, Eq. (42) gives t 
be rewritten as 



R 
23 



■1/t + e and t ^ {t + l/t)xc§3- The latter can 



ms M- 



0(1) 



(43) 



m2 V 

Notice that in this case cannot be arbitrarily small, otherwise M3 would exceed 



the GUT or Planck scale. In addition 

sex 



R 



sx Xc -1, 



(44) 



The large masses m2 and mi are dominated by A^i and N2 with competitive strenght^. 
In particular, the first row of R determines leptogenesis, as we are going to discuss. 
The smallest mass ma is instead dominated by A^3, which decouples the more x is 
small. The Yukawas are 



SO that 



M3 



mi 



u2i 



M2 



-a¥* + o(x6)^^u:, , 

(45) 



' + ,/^U*UiWtW3 ) , 



^4 = Tr^\'^\h^ Im( -JxUiWtlWsl^ + JxUsW^r 
iV2@ M@ y mi \ 

which turns out to be proportional to Uez- 

If c{| = e < 1, one has t2^ — t{l — et) and t ^ et (i + 1/i )xc^, which in turn implies 



(46) 



ms / M3 



1712 V M@ 



0{e) 



(47) 



In this case m3 is allowed to be small and the limit m3 — > can be applied. In 
addition 



R^ 



s c 

-c s ^ -X 1 +C»(e^ex,x')■ 
cx -Xs -1 



(48) 



The Yukawas are 

^i/3i = 



M3 



V "^1 



M2 



-a¥:-o(x)J— u*) , 

V mi 



(49) 



^Assuming for instance ai = and tri-bimaximal mixing one has c = -\/2/3, s = -\/l/3; while for 
Q!2 — Q!i = TT one has c = v^, s = —i. 
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so that 



M2 M3 ,^2 /^T.__ f ~ TTT T,Tt7*lT,-!,r 12 , - TTT ^-.7*1^^7 |2 , /^l 



(50) 

which is again proportional to Uez and suppressed by a factor of e with respect to 
(46). Hence, it vanishes in the hmit ma = due to the fact that the second and 
third rows of are ahgned. The contributions from /J^^^'^^^^ are neghgible even for 
the future sensitivities. 



Summarizing, in the inverted-hierarchical case can be at hand of future experiments 
only if M2 ~ M@, m'i/m@ ~ Mq/M^ and \Ue3\ is large. The dependence on 02 and S is 
quite complicated and also ai plays a role. 

With an IH light neutrino spectrum, the RGE effects might be important, in particular 
for the solar angle and mass squared difference. Nevertheless, the above conclusions 
remain valid. In practice, one has just to rescale m and ttiq at high energy by a numerical 
factor, and check whether the parameters of U, when evolved at low energy, fall within 
the experimental window. It is well known [39] that this is always the case if |q;2 — q;i| — tt, 
namely if the solar pair of eigenstates are of pseudo-Dirac type. Then, /J^ and all the 
Cij are mildly dependent on a2 and essentially depend only on S. This will turn out to 
be the case also for Yg. 

Wc illustrate this considering an example of the first case where the RGE effects turn 
out to be quite small. Wc consider tan /5 = 30 and take at Mpi: rh = 1.5(1, 0.97, 10^^)m@, 
(M3,M2) = (8,0.1)M@, Ml = 10" GeV, 623 = 45°, ^12 = 35°, ^13 = 10°, 6 = 37r/2, and 
any q;2 = ai + tt. At low energy, the neutrino spectrum is viable, tti = (1, 0.98, 10~^)m@, 
and the angles of the MNS are 623 = 44.4°, 9i2 = 35°, ^13 = 10.4°. We obtain /f^ = 
-4.8, = -0.45, C31 = 1.4, C32/4 = C21 = 3 X 10"^. Hence, for the point P, 
= _5 X 10-^rff With tan p = 30, the FV contribution to 4 is at hand of future 
experiments = O.Oldg^^; as for the BR of r — 67, r fi'-f, fi ej, they are smaller 
than their corresponding experimental limits by factors of 0.7, 4 x 10~^, 0.15, respectively. 

Fig. 4 shows the dependence of and Yb on d, keeping for the remaining parameters 
the same set as before. The curves are absent for 6 very small or close to 2n, where 
C31 > 10 because some of the Yukawas in Y^ blow up. Since both Yb and Jj^ are directly 
proportional to \Ue3\, the plot can be adapted correspondingly to other values of 6*13. 
For a sizeable de and Yb, ^13 cannot be smaller than a few degrees. The value chosen 
for Ml is 10^^ GeV which corresponds to the regime of temperatures where all lepton 
flavours, but the electron one, are in thermal equilibrium, since tan (3 — 30. In this case, 
Eq. (14) applies. We have checked that for smaller values of tan /3, where only the tau is 
in equilibrium, the final baryon asymmetry does not change significantly. 

In Fig. 5 we show the dependence of Jj^ and Yb on 7713/7712 (left and right panels, re- 
spectively), while keeping m3/m2 — 0.08 M@/M3 as in the previous example and assuming 
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Figure 4: Dependence of Jf^ and Yb on 6 for 02 = + tt: we take 0:2 = 0, tt/2, tt, 37r/2. 
The vertical light-grey regions are excluded by requiring perturbative Yi,. For the choice 
of the other parameters, see the text. 




mi/rrih 



Figure 5: Dependence of Jj^ and on rrii/mY^. For normal hierarchy mi/mh = mi/m^, 
while inverted hierarchy mi/mh = m^/m2- For the choice of the other parameters, see 
the text. The dashed lines are the result of the 1-flavour approximation. 

the same set for the remaining parameters (in particular we selected 0:2 = 0, ai = tt). On 
the right-hand side we compare the baryon asymmetry computed with the flavour effects 
included to the baryon asymmetry computed in the one-single flavour approximation. We 
see that there is almost one order of magnitude difference, both in the normal and in the 
inverted hierarchical case. The main reason is that in the one-flavour case the total lepton 
asymmetry is strongly washed out since the wash-out parameter m = rhf. + + m,- is 
larger than m@, yielding a suppressed baryon asymmetry. On the contrary, in the flavour 
approximation, the asymmetries in the electron and muon flavours are only weakly washed 
out. Furthermore, in the inverted hierarchy case, the one- flavour approximation leads to 
a suppression in the CP asymmetry that goes like mQ/m@, while including flavours the 
individual CP asymmetries go as m@. 



18 



5 Conclusions 



The predictions for the electron EDM in the context of the supersymmetric seesaw and 
mSUGRA have been investigated. First, we showed the existence of an indirect upper 
bound on de from negative searches for r — > 67, within the supersymmetric seesaw. This 
indirect bound may be even stronger than the present experimental direct upper limit. 
Planned searches for rfg, improving the sensitivity by about three orders of magnitude, 
will be able to supersede the indirect bound from r ^ 67 and provide considerable tests 
of the seesaw-induced effects. We identified in a model-independent way the patterns of 
seesaw models that lead to a potentially observable electron EDM, considering in turn the 
case of normal and inverted light neutrino spectra. A widely sphtted spectrum for right- 
handed neutrinos is a crucial ingredient, as well as the relation mi/mh ~ M@/Ms, where 
M@ = v^/m®. Indeed, in the limit of a vanishing lightest neutrino mass, the electron 
EDM drops much below the planned sensitivities. 

The seesaw interactions may also be responsible for the generation of the baryon 
asymmetry of the Universe via the mechanism of leptogenesis. The patterns of seesaw 
models identified requiring an electron EDM within future experimental sensitivities allow 
to extract the value of Mi required to generate a sufficiently large baryon asymmetry via 
thermal leptogenesis. The importance of taking into account fiavour effects has been 
emphasized. Our findings show that a large enough baryon asymmetry may be achieved 
through thermal leptogenesis for those patterns which give rise to a large electric EDM 
and suppressed r — > 67, r — > /X7, — > 67 branching ratios. However, a sufficiently 
large baryon asymmetry is reached only for large values of Mi , at least larger than about 
lO^-*^ GeV. Since we are dealing with a supersymmetric leptogenesis set-up, we should 
face the problem arising from the so-called gravitino bound. The latter is posed by the 
possible overproduction of gravitinos during the reheating stage after infiation, see for 
instance [40]. Being only gravitationally coupled to the SM particles, gravitinos may 
decay very late jeopardising the successfuU predictions of Big Bang nucleosynthesis. This 
does not happen, however, if gravitinos are not efficiently generated during reheating, 
that is if the reheating temperature Trh is bounded from above, Trh < 10^° GeV [40] . 
The severe bound on the reheating temperature makes the generation of the RH neutrinos 
problematic (for complete studies see [19, 32]), if the latter are a few times heavier than 
the reheating temperature, rendering the thermal leptogenesis scenario unviable. 

In view of the above, if a large de is measured in the near future and assuming that 
its main contribution comes from CP-violating effects induced by the see-saw Yukawas, 
either the baryon asymmetry is not explained within the thermal leptogenesis scenario 
or leptogenesis occurs in a non-thermal way. This second alternative stems from the 
fact that the RH neutrinos might be generated not through thermal scatterings, but by 
other mechanisms, for example during the preheating stage [41], from the inflaton de- 
cays [41, 42] or quantum fluctuations [43]. In these cases, the baryon asymmetry depends 
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crucially on the abundance of RH neutrinos and sneutrinos generated non-thermally. For 
instance, these heavy states may be produced very efficiently during the ffist oscillations 
of the inflaton field during preheating up to masses of order (10^^ — 10^^) GeV [41]; the 
final baryon asymmetry is generated of the right-order of magnitude if the flavour lepton 
asymmetries satisfy the mild condition > 10-^(10^^ GeV /Trh) {Mi/ 10^^ GeV), as one 
can readily deduce from Ref. [41]. 
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A Constraints from LFV e EDM 



In Fig. 6 we display Q/, If^^ and I^^^^ from the present 90% C.L. limits [44], BR(/i 
67) < 1.2 X 10-^^ , BR(r fx-f) < 6.8 x 10"^ , BR(r 67) < 1.1 x 10"^. The upper 
bounds Jg^^*-^ and /"^^^ are also shown taking de < l-Qx 10^^'^ e cm. We consider mSugra 
with oo = rriQ + M1/2. For better sensitivities, the values of Cfj' have to be multiplied by 

a factor ^BR^^VBRP^ while those of Jf and Jf by a factor 4"VrfP^ We recall that 

Ml ~ O.4M1/2 and m|, ~ + 0.15M^y2- The plots are adapted from those in [14], where 
the reader can find more details and references to the literature. 
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Figure 6: Contours of C^j (upper plots), from the present 90% C.L. limits [44]. The lower 
plots show the contours of I^^^'^ and J"*^^^ taking < 1.6 x 10~^^ e cm. 
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